Accurate simulations of planar topological codes cannot use cyclic boundaries 
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Cyclic boundaries are used in many branches of physics and mathematics, typically to assist 
the approximation of a large space. We show that when determining the performance of planar, 
fault-tolerant, topological quantum error correction, using cyclic boundaries leads to a significant 
underestimate of the logical error rate. We present cyclic and non-cyclic surface code simulations 
exhibiting this discrepancy, and analytic formulae precisely reproducing and explaining the observed 
behavior in the limit of low physical error. 



Topological quantum error correction (TQEC) is the 
lowest overhead technique for achieving reliable large- 
scale quantum computation given a 2-D lattice of qubits 
with nearest neighbor interactions pQ. As such, accu- 
rate simulations of the performance of such codes under 
physically realistic assumptions are of great practical rel- 
evance. Much of the existing literature on TQEC uses 
cyclic boundaries (2H3]- While this is sufficient to obtain 
the threshold error rate pth of a given TQEC scheme, 
we show that this is not sufficient to obtain an accurate 
logical error rate at a given code distance d and physical 
error rate p < p t h- Indeed, we show that cyclic bound- 
aries lead to an exponentially growing underestimate of 
the logical error rate with d. 

Without loss of generality, we shall focus on the sur- 
face code [rJHTU] making use of minimum weight perfect 
matching for decoding |llH14j . Our results are relevant 
to any periodic, fault-tolerant implementation of any pla- 
nar TQEC code. An up-to-date and detailed review of 
the surface code can be found in [TS]. We shall use the 
error models (standard depolarizing) and circuits (8-step 
error detection cycle using only Z basis initialization and 
measurement) described in this review. 

The patterns of stabilizers [TH] associated with distance 
d = 3 non-cyclic and cyclic surface codes are shown 
in Fig. [TJ Consider Fig. [2] This shows how errors 
propagate and create detection events both away from 
and near boundaries of the qubit array (the latter only 
when assuming non-cyclic boundaries) . It is conceptually 
straightforward to analyze the propagation of all possible 
single errors and determine the total probability of any 
given pair of detection events due to all distinct single er- 
rors. We can visualize each such probability as a cylinder 
with space-time cndpoints corresponding to the detection 
events and diameter proportional to the total probabil- 
ity. We call such a visualization a nest, and each cylinder 
a stick. Nests for X stabilizer measurement in the sur- 
face code assuming non-cyclic and cyclic boundaries are 
shown in Figs. [3] and |4j respectively. These figures were 
generated using our tool Autotune |17j . modified to per- 
form a simple sum of probabilities rather than the more 
complex polynomial used in prior work |18) . Due to the 
complexity of these figures, their raw data has been in- 
cluded in the Supplementary Material. 

The probability of each type of logical error per round 
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FIG. 1: Patterns of stabilizers of a distance d = 3 surface 
code with a. non-cyclic boundaries, b. cyclic boundaries. Ex- 
amples of logical operators of each type have been given. 



of error detection for various distances d and non-cyclic 
and cyclic boundaries as a function of the physical depo- 
larizing error rate p is shown in Figs. [5} |10| It can be seen 
that the logical error rates of the cyclic case are signifi- 
cantly lower, particularly logical X\ and Z 2 errors. Wc 
can determine the low p asymptotic forms of the curves 
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FIG. 2: 2-D surface code (grey). Time runs vertically. 
Squares represent initialization to jO), circles represent Z basis 
measurement. Slashed squares represent initialization to +), 
slashed circles represent X basis measurement (both achieved 
using Hadamard gates), a. A single error leading to a pair 
of detection events (green ellipses). A detection event is a se- 
quential pair of measurements with differing value. Red lines 
show the paths of error propagation, b. An error leading to a 
single detection event due to proximity to a boundary of the 
lattice. 





FIG. 4: Nest of 7 rounds of Fig.JTJs fault-tolerant X stabilizer 
measurements. 



d 


A x 


A z 


A Xl 


A Zl 


4 


3.97 x 10 2 


4.70 x 10 2 


1.23 x 10 2 


2.76 x 10 2 


6 


1.67 x 10 4 


2.09 x 10 4 


1.97 x 10 3 


6.64 x 10 3 


8 


7.02 x 10 5 


9.34 x 10 5 


2.94 x 10 4 


1.49 x 10 5 


10 


2.93 x 10 7 


4.18 x 10 7 


4.23 x 10 5 


3.21 x 10 6 



TABLE I: Low p analytic asymptotic formulae derived from 
first principles, not simulation. Each entry represents a curve 
of the form Ap d ' 2 . The definitions of the logical operators X, 
Z, Xi, Z\ can be found in Fig. [I] 




FIG. 3: Nest of 7 rounds of Fig.JT^i fault-tolerant X stabilizer 
measurements. 



in Figs. [5] fusing the work of 191 . These analytic expres- 
sions are shown as dashed lines and detailed in Table H] 
and were derived completely independently of the simula- 
tions. The agreement with the simulation data is perfect 
to numerical precision, implying the processes leading to 
failure in the non-cyclic surface code are well understood. 
As the non-cyclic code distance d increases, the num- 



ber of minimum stick topologically nontrivial paths con- 
necting opposing boundaries grows exponentially (con- 
sider larger versions of Fig. [3]). Such paths can be associ- 
ated with logical operators of minimum weight d, each of 
which crudely speaking increases the probability of log- 
ical error. More precise discussion can be found in |19j . 
By contrast, in a periodic, fault-tolerant, distance d sur- 
face code with cyclic boundaries, there are precisely d 
minimum stick topologically nontrivial cyclic paths. Di- 
agonal sticks never zigzag and hence never create addi- 
tional cyclic paths. This is a generic property of periodic, 
fault-tolerant TQEC. Logical errors will occur 50% of the 
time when physical errors occur corresponding to half the 
sticks in a given minimum weight logical operator. This 
implies the low p asymptotic logical error rate, for even 
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FIG. 5: Probability of logical X error as defined in Fig.[T^, as 
a function of the depolarizing error rate p for various distances 
d. The analytic forms of the dashed curves can be found in 
Table |U 



FIG. 7: Probability of logical Xi error as defined in Fig.JTJj as 
a function of the depolarizing error rate p for various distances 
d. The analytic forms of the dashed curves can be found in 
Table |U 
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FIG. 6: Probability of logical Z error as defined in Fig.[]J, as a 
function of the depolarizing error rate p for various distances 
d. The analytic forms of the dashed curves can be found in 
Table H] 



FIG. 8: Probability of logical Z\ error as defined in Fig.fTJa as 
a function of the depolarizing error rate p for various distances 
d. The analytic forms of the dashed curves can be found in 
Tabled] 



distances, will be 




where e is the probability of a stick within a minimum 
stick logical operator. Note that the symmetry of a cyclic 
boundary surface code implies that all such sticks will 
have the same probability. 

From the Supplementary Material, we find that a Z\ 
logical operator (left to right) is associated with sticks 
of probability e = 4.8p, for example by considering the 



(4, 1, 6) to (4, 3, 6) stick. A Z 2 logical operator (top to 
bottom) is associated with sticks of probability e = 3.2p. 
Inserting these expressions into Eq. [I] and noting that 
by symmetry the logical X% and Z2 error rates will be 
the same, as will the logical X 2 and Z\ error rates, we 
can see in Table [I] that at distance d = 10 the cyclic 
logical Z error rate is already over a factor of 10 too 
low and the cyclic logical X error rate nearly 2 orders of 
magnitude too low. The degree of underestimation grows 
exponentially with d. 

In summary, we have shown that one of the most 
commonly used simplifying assumptions in physics and 
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mathematics, namely cyclic boundaries, cannot be used 
if one wishes to accurately study the logical error rates 
of planar topological codes. The level of inaccuracy in- 
troduced by this seemingly harmless assumption is expo- 
nentially large in the code distance d. This implies that 
cyclic boundaries should be avoided when performing ac- 
curate studies of practical topological quantum error cor- 
rection. We hope that this work encourages the quan- 
tum information community to seriously question the ac- 
curacy of results derived using unphysical assumptions. 
Two common assumptions that we feel deserve partic- 
ular scrutiny are assuming perfect multi-body quantum 
measurements, and assuming arbitrarily long range co- 
herent interactions without time, overhead, or error rate 
penalty. 



FIG. 9: Probability of logical X2 error as denned in Fig.[T]D as 
a function of the depolarizing error rate p for various distances 
d. The analytic forms of the dashed curves are identical to 
those of Fig. [8] and can be found in Table [I] 




Tx10" 5 1x10" 4 1x10" 3 1x10" 2 

Depolarizing probability (p) 



FIG. 10: Probability of logical Z2 error as defined in Fig.fTJs as 
a function of the depolarizing error rate p for various distances 
d. The analytic forms of the dashed curves are identical to 
those of Fig. [7] and can be found in Table [I] 
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Appendix A: Supplementary Material 

We include below the raw data for the nests of Fig. [3] 
and Fig. [4] The initial two groups of three numbers in 
each line are the space-time coordinates of a cylinder, 
with the i, j coordinates corresponding to the numbers 
m Fig.[TJ and the* coordinate corresponding to the num- 
ber of rounds of error detection since the first, doubled. 
The time t therefore takes the values {0, 2, ... , 12}. The 
doubling of the time coordinate is purely to improve the 
symmetry of the nests. The final number in each line is 
the diameter of each cylinder and is the total probabil- 
ity of detection events at the end points of each cylinder 
from single errors when using a standard balanced depo- 
larizing error model of strength p = 0.04. This value of p 
was chosen simply to produce a clear nest. The diameter 
of each cylinder scales linearly with p. 

1. Data for Fig. [3] 
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